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VALUES  OF  THE  PHOTOMETRIC  PARAMETERS  OF  MARS 
AND  THEIR  INTERPRETATION 

By  W.  R.  Weaver,  W.  E.  Meador,  and  G.  P.  Wood 
Langley  Research  Center 

SUMMARY 

Photometric  data  for  the  Martian  disk  obtained  by  Thorpe  from  television- 

V 

camera  pictures  taken  by  Mariner  9  have  been  used  to  derive  values  for  the 
parameters  in  the  photometric  function  due  tw>  Meador  and  Weaver.  The  photo¬ 
metric  function  can  thus  be  applied  to  Mars  and  to  the  design  of  cameras  for 
photography  of  Mars.  Values  of  the  parameters  are  for  one  wavelength  only, 

0.56  pm,  and  are  average  or  effective  values  for  the  Martian  disk.  The  values 
derived  herein  have  been  interpreted  to  provide  semi -quantitative  information 
on  the  soil.  The  results  do  not  disagree  with  other  indications  that  the 
effective  mean  diameter  of  the  surface  particles  is  about  400  pm.  Another  re¬ 
sult  is  that  the  mean  intercenter  spacing  of  adjacent  particles  may  be  about  4/3 
of  the  mean  diameter. 


INTRODUCTION 

It  would  be  desirable  to  have  available  a  photometric  function  for  plane¬ 
tary  surfaces  that  would  describe  the  brightness  of  a  surface  in  reflected  sun¬ 
light  as  function  not  only  of  the  angles  of  incidence,  emission,  and  phase,  but 
also  as  function  of  the  physical  properties  of  the  soil  that  forms  the  surface. 
The  commonly-used,  indeed  universally-used,  photometric  ftinction  for  correlating 
and  analyzing  data  on  planetary  surfaces  and  for  describing  the  reflection  from 
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(i  r  brightness  of)  planetary  surfaces  is  the  well-known  Minnaert  function. 

T!  is  function  attempts  to  describe  the  brightness  as  it  depends  on  the  geometry 
of  the  situation;  namely,  the  angles.  The  Minnaert  function,  however,  has  very 
insufficient  basis  in  theory,  is  not  derived  from  first  principles,  is  restrict¬ 
ed  in  validity  to  certain  classes  of  scattering  geometries,  and,  of  consider¬ 
able  importance,  can  not  be  used  to  interpret  photometric  behavior  in  terms  of 
the  physical  properties  of  the  reflecting  surface. 

This  undesirable  situation  has  now  been  improved.  W.  E.  Meador  and  W.  R. 
Weaver  of  the  Langley  Research  Center  have  derived  a  new  photometric  function 
in  an  as  yet  unpublished  paper  entitled  "A  Photometric  Function  for  Diffuse 
Reflection  by  Particulate  Materials.”  This  photometric  function  is  not  only 
grounded  in  theory  but  has  been  verified  for  several  soils  in  laboratory  experi¬ 
ments.  It  is  applicable  to  all  geometries  of  scattering,  it  takes  into  account 
single  scattering,  multiple  scattering,  and  interparticle  shadowing,  and  it 
also  contains  parameters  the  values  of  which  depend  on  the  physical  properties 
of  the  surface.  The  investigation  on  which  this  new  photometric  function  is 
based  is  not  yet  complete,  inasmuch  as  interpretation  of  all  but  one  of  these 
parameters  in  terms  of  physical  properties  of  the  soil  is  not  yet  on  a  very 
exact  basis,  but  the  investigation  is  at  a  stage  that  permits  qualitative 
and  semi-quantitative  statements  about  the  size  and  the  albedo  of  the 
particles  that  make  up  the  soil  and  the  degree  of  packing  of  the  particles,  or 
the  porosity  of  the  soil. 

The  purpose  of  this  paper  is  to  use  this  new  photometric  function  together 
with  the  latest  and  apparently  the  best  data  on  the  photometric  properties  of 
Mars  to  obtain  both  a  complete  statement  of  the  whole-disk  photometry  of  the 
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planet  (at  one  wavelength)  and  some  at  least  qualitative  information  on  the 
effective  average  properties  of  the  soil'. 

THE  PHOTOMETRIC  FUNCTION 

The  photometric  function  derived  by  Meador  and  Weaver  is 

=  TuySpiffri  .  cos  e)  [(Ua0  cos  a)  f(i'e'“;a2) 

+  a^  (cos  i  +  cos  e)J  (1) 

where  $(i,e,a)  is  the  brightness  (normalized  to  unity  at  i  *  e  =  0) ;  i  is 
the  angle  of  incidence;  e  is  the  angle  of  emission;  a  is  the  phase  angle;  and 

v 

aQ,  a^,  are  empirical  parameters  that  contain  information  about  the  surface. 
The  factor  f(i,e,a;a2)  is  given  by 


f(i,e,a;a2)  -  e 


v/1  expii 
0  < 
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y  -  |3ir(2g-l)x  +  6xsin  1x  +  2(2+x2)  (1-x2)  dx 


(2) 


where 


4a2  (1  ♦  cos  a) 
y  3  sin  a 

ira2  (cos  i  +  cos  e)  [cos  Y  +  cos  (o-Y))  cos 

v  =  - 

sin  a  cos  i  cos  e 

cos  a  =  cos  i  cos  e  +  sin  i  sin  e  cos  $ 

(cos  i  »  cos  e)  cos  Y  cos  (a-Y)  cos  ^ 
g  [cos  Y  +  cos  (a-Y)]  cos  i  cos  e 

and  ip  is  the  an^isJ  between  the  surface  normal  and  the  scattering  plane,  <J>  is 
the  azimuthal  angle  between  the  plane  of  incidence  and  the  plane  of  emission, 
and  Y  is  the  angle  between  the  incidence  direction  and  the  projection  of  the 
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an 
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surface  normal  on  the  scattering  plane.  The  complete  normalized  photometric 
function  is  the  combination  of  equation  (1)  with  the  f-factor  of  equation  (2). 

THE  DATA 

In  reference  1,  Thorpe  reported  the  results  of  a  statistical  comparison  of 
more  than  2500  Mariner  9  television-camera  pictures  of  Mars  that  was  made  to 
determine  the  integral  (as  opposed  to  regional)  photometric  properties  of  the 
planet.  He  reported  his  results  in  the  form  of  the  Minnaert  parameters  BQ 
and  k.  At  this  point  we  remind  the  reader  that  the  Minnaert  function  is 

$  =  Bn  (cos  i)k(cas  e)k  1 

and  that  a  Minnaert  plot  is  a  plot  of  Jtn(4>  cos  e)  against  K,n(cos  i  cos  e)  . 
Thorpe  obtained  the  values  of  the  Minnaert  parameters  for  a  portion  of  the 
mission  that  followed  abatement  of  the  dust  storm  (revolutions  of  Mariner  9 
nos.  220  to  260).  The  values  of  these  parameters  are  shown  in  Thorpe's  Table  I. 
The  Minnaert  coefficient  Bq  has  been  plotted  in  the  present  paper's  figure  1. 
Also  shown  in  figure  1  is  the  BQ  obtained  from  measurements  in  the  laboratory 
by  Meador  and  Weaver  on  Colorado  basalt  (mafic  latite  porphyry)  of  size  150  to 
300  ym.  The  values  of  Bq  for  both  curves  are  relative,  as  there  is  no  means 
for  normalizing  them  to  the  same  standard  other  than  by  simple  matching  at  some 
point.  The  curve  for  Colorado  basalt  was  adjusted  to  agree  with  the  Mariner  9 
curve  at  a  phase  angle  of  50°.  These  curves  indicate  that  the  trends  in  the 
plot  of  Thorpe's  data  are  in  reasonable  agreement  with  laboratory  results. 

Thorpe's  values  for  the  Minnaert  exponent  k  have  been  plotted  in  figure 
2.  Also  in  figure  2  are  plots  of  k  for  two  sizes  of  Colorado  basalt  as  deter¬ 
mined  by  Meador  and  Weaver  in  their  laboratory  investigation.  They  are  shown  in 
order  to  indicate  whether  the  magnitude  and  the  trend  of  Thorpe’s  k  appear  to 
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bt  in  reasonable  agreement  with  laboratory  results.  They  do.  Also  shown  in 
figure  2  are  values  of  k  from  Binder  and  Jones  (reference  2)  (corrected  from 
an  erroneous  plot)  for  regions  at  and  near  Chryse,  Tharsis,  and  Amazonis,  and 
values  from  Young  and  Collins  (reference  3)  for  Ophir,  Aeolis,  and  Solis  Lacus. 

DERIVATION  OF  PARAMETERS 


The  theoretical  expression  of  Meador  and  Weaver  for  $  and  Thorpe's  photo¬ 
metric  data  on  Mars  were  used  together  to  derive  numerical  values  for  the  para¬ 
meters  a^,  a^,  and  a2  in  the  Meaddr-Weaver  photometric  function.  An  itera¬ 
tive  procedure  was  followed  that  makes  use  of  the  fact,  derived  by  Meador  and 

% 

Weaver,  that  the  Minnaert  function  is  a  good  approximation  to  the  Meador-Weaver 
function  for  three  geometries.  One  of  these  is  a  geometry  for  which  the  angular 
deviation  from  the  mirror  point  is  sufficiently  small  that  the  corresponding 
plot  is  linear.  The  constant  part  of  the  slope  of  a  Minnaert  plot  of  equation 
(1),  which  is  k,  is  calculated  from  the  equation 


k(an,.,,a„cO  =■  j*  C°S  e)lsl| 

012  £+q  |d[£n  (cos  i  cos  e)]jo  = 


(3) 


constant 

where  £  =  angular  deviation  from  mirror  geometry.  The  iterative  procedure  is: 

1.  A  value  for  a^  between  -1  and  +1  is  selected. 

2.  Values  of  a^  and  a2  are  determined  such  that  equation  (3)  agrees 
with  Thorpe's  data  for  k  as  a  function  of  a. 

3.  Thorpe's  values  for  brightness  are  then  normalized  to  match  $  for 
a  =  50°  and  e  =  0°. 

4.  Thorpe's  values  of  normalized  brightness  are  graphically  compared  with 
the  4>  of  equation  (1)  for  the  scattering  geometries  e  =  0°  and  varying  a. 
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(brightnesses  corresponding  to  these  geometries  have  been  shown  to  be  well 
approximated  by  Minnaert  functions.) 

5.  Then  other  values  of  a^  are  chosen  and  the  entire  procedure  is 
repeated  until  a  plot  of  $  vs  a  agrees  closely  with  a  plot  of  Thorpe's 
normalized  brightness  vs  a. 

The  results  of  this  matching  of  theory  and  experiment  are  that 

ag  =  0.55 
a^  =  0.60 
&2  =  0.18 

To  see  that  equation  (1),  when  used  with  these  values  for  the  parameters,  pro¬ 
vides  good  agreement  with  Thorpe's  data,  figure  3  may  be  examined.  The  of 
equation  (1)  is  shown  by  the  curve,  Thorpe's  data  by  the  circles.  Figure  4 
shows  the  agreement  between  the  effective  Minnaert  exponent  from  equation  (3) 
and  the  values  of  Thorpe. 

The  above  values  of  the  parameters  a^  together  with  equation  (1)  permit 
calculations  of  the  brightness  of  Mars  for  any  geometry.  They  also  provide  in¬ 
formation  for  use  in  designing  cameras  for  photography  of  Mars. 

INTERPRETATION  OF  VALUES  OF  PARAMETERS 

The  development  of  the  photometric  function  of  Meador  and  Weaver  has  lead 
to  the  ability  to  make  qualitative,  and  even  approximate  quantitative,  interpre 
tations  of  the  numerical  values  of  the  three  parameters  in  the  photometric  func 
tion  in  terms  of  the  physical  properties  of  the  surface.  Because,  however, 
further  theoretical  and  experimental  investigation  of  the  exact  meaning  of  each 
of  these  parameters  is  needed,  and  because  of  the  inadequacy  of  data  about  Mars 
deductions  presented  herein  about  the  physical  properties  of  the  surface  mater- 
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iul  should  be  considered  to  be  tentative  and  to  be  relative  rather  than  abso¬ 
lute  indications  of  properties  of  the  surface  material  of  Mars. 

The  parameter  ag  is  permitted  to  have  any  value  between  the  limits  of 
-1  and  +1.  An  excess  of  forward  scattering,  as  could  occur  by  means  of  trans¬ 
mission  through  small  particles,  is  indicated  by  negative  values  of  ag.  An 
excess  of  backscattering  is  indicated  by  positive  values  of  ag.  Isotropic 
scattering  is  indicated  by  ag  equal  to  zero.  The  determination  of  ag  from 
experimental  data  should  yield  some  concept  of  mean  particle  size  through  the 
relation  of  that  property  to  the  particle  transparency. 

In  interpreting  the  values  of  the  parameters  for  Mars,  it  is  helpful  to 
compare  the  values  with  those  measured  for  other  materials  in  the  laboratory. 

In  Table  I  are  shown  the  results  of  some  measurements  by  Meador  and  Weaver  on 
two  size  ranges  of  each  of  two  different  materials.  Notice  that  as  the  size  of 
the  laboratory  materials  increases  the  value  of  aQ  increases.  The  value  of 
aQ  for  Mars  of  0.55  indicates  that  the  particles  are  sufficiently  large  that 
backscattering  predominates  over  forward  scattering.  If  the  material  on  the 
surface  of  Mars  has  values  of  permittivity,  permeability,  and  electrical  conduc¬ 
tivity  that  are  not  too  different  from  those  of  the  laboratory  materials,  then 
the  large  positive  value  of  ag  implies  a  mean  particle  size  considerably  larger 
than  225  ym. 

The  parameter  aj  is  proportional  to  the  ratio  of  multiple  to  single 
scattering  and  thus  relates  to  the  single-particle  albedo.  The  limits  of  a^, 
are  0  and  °°,  corresponding,  respectively,  to  the  limit  of  very  dark,  single¬ 
scattering  materials  and  to  the  limit  of  very  bright,  Lambert  surfaces.  Gener¬ 
ally,  as  particle  size  is  decreased  some  of  the  radiation  that  would  be  absorbed 
by  the  larger  particles  is  transmitted  through  the  smaller  ones  to  be  scattered 
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a^ain.  This  process  enhances  multiple  scattering  and  is  the  explanation  for 
the  increase  in  with  decrease  in  particle  size  shown  in  Table  I.  Since 

the  transmitted  and  rescattered  light  partially  emerges  as  reflected  light, 
this  process  also  accounts  for  the  common  observation  that,  for  sufficiently 
small  particles,  surfaces  are  brighter  for  smaller  particles  of  the  same  compo¬ 
sition  than  they  are  for  larger.  The  value  of  a^  of  0.6  indicates  only  that 
multiple  scattering  is  a  significant  contributor  to  the  brightness.  Thus  any 
light-scattering  model  of  the  Martian  surface  must  take  into  account  both  single 
and  multiple  scattering  to  be  adequate. 

3 

The  parameter  a„  is  a  packing  factor.  It  has  the  value  n(cd)  where 
d  is  the  mean  particle  diameter  and  n  is  the  particle  number  density.  To  the 
extent  that  radiation  passing  through  the  gaps  between  particles  is  undiffracted, 
the  value  of  c  lies  between  0.5  and  1.  We  choose  for  use  here  the  value  0.75. 
Thus,  with  a2  ■  0.18  for  Mars,  the  mean  intercenter  particle  distance  between 
nearest  neighbors,  which  is  the  cube  root  of  1/n,  is  indicated  to  be  4/3  of  the 
mean  diameter  of  the  particles.  This  result  should  be  regarded  as  a  lower  limit 
because  of  the  expectation  that  c  is  an  increasing  function  of  wavelength 
throughout  the  visible  and  infrared  range. 

The  lower  limit  of  a 2  approaches  zero  and  the  upper  limit  approaches 
unity.  The  smaller  the  value  of  a2  the  more  porous  is  the  structure  of  the 
surface.  That  the  value  of  a2  for  Mars  of  0.18  is  approximately  midway  in 
the  range  of  a2  for  the  laboratory  materials  indicates  that  the  surface  of 
Mars  is  packed  to  about  the  same  extent  as  the  surface  of  the  laboratory  mater¬ 
ials.  The  laboratory  materials  were  applied  by  sifting,  or  airdropping  the 
particles  to  the  surface  from  a  height  of  about  0.5  m.  The  degree  of  packing 
that  occurred  is  not  inconsistent  with  a  Martian  surface  formed  by  strong 
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a>  olian  forces.  A  cautionary  note,  however,  is  needed.  The  measure  of  packing 
derived  herein  applies  only  to  that  layer  of  surface  material  that  contributes 
to  photometry—  probably  only  a  few  particle  diameters.  The  use  of  radiation  of 
much  longer  wavelength  than  that  of  visible  light;  e.g.,  radar  wavelengths,  is 
needed  to  determine  the  bearing  strength  of  the  soil. 

OTHER  INDICATIONS  OF  PARTICLE  SIZE 

Reference  4  reports  the  results  of  an  investigation  of  the  thermal  proper¬ 
ties  of  portions  of  the  surface  of  Mars  by  an  infrared  radiometer  that  was 

carried  on  Mariner  9.  Figure  5  of  this  reference  shows  that  as  the  effects  of 

♦ 

the  dust  storm  decreased,  the  best-fit  planet-wide  value  of  the  thermal  inertia 

-2  -1  -1/2 

approached  the  value  0.007  ±.  0.001  cal  cm  K  sec  .  Figure  11  of  the  same 

reference  shows  that,  on  the  basis  of  the  observations  made  with  the  radiometer 
at  a  latitude  of  -30°  and  with  a  reasonable  assumption  about  the  value  of  the 
volume  specific  heat,  the  above  value  of  thermal  inertia  indicates  a  mean  parti¬ 
cle  diameter  of  360  um.  The  statement  is  also  made  that  these  results 

are  not  strong  functions  of  latitude,  date,  or  the  predawn  time  of  day  at  which 
the  surface  temperatures  were  measured. 

There  is  one  further  piece  of  evidence  about  the  size  of  the  surface  parti¬ 
cles  on  Mars.  As  is  well  known,  the  minimum  value  of  the  speed  of  the  wind  on 
Mars,  at  the  top  of  the  boundary  layer,  that  is  required  to  initiate  the  motion 
of  surface  particles  is  rather  high,  about  50  m/s  (ref.  5).  Instances  where 
such  motion  occurs  seem  to  be  rather  frequent.  One  would  expect,  therefore, 
that  the  initial  motion  of  particles  would  be  that  of  particles  for  which  the 
minimum  value  of  wind  speed  is  sufficient  and  that  a  significant  portion  of  the 
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soil  is  composed  of  such  particles.  As  is  shown  in  reference  5,  the  particle 
diameter  for  which  the  required  wind  speed  is  a  minimum  is  about  450  pm. 

There  is  thus  no  disagreement  between  other  indications  of  particle  size 
and  the  implication  from  the  value  of  a^  that  the  size  should  be  considerably 
greater  than  225  ym. 
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TABLE  I.-  Photometric  Parameters  from  Laboratory  Measurements 


Material 


Colorado  basalt 


Basalt  dune  sand 


Size 

range,  ym 

Mean 
size,  ym 

fo 

60  -  150 

105 

-0.4 

0.28 

0.32 

150  -  300 

225 

-0.1 

0.26 

0.15 

50  -  200 

125 

-0.1 

0.25 

0.17 

100  -  320 

210 

+0.05 

0.20 

0.09 

Figure  L  -  Comparison  of  the  Minnaert  coefficient  for  Mars  (Mariner  9,  Thorpe) 
and  a  basalt 
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Figure  2.  -  Comparison  of  the  Minnaert  exponent  for  Mars  (derived  from  Earth-based 
and  Mariner  6,  7,  and  9  data)  and  a  basalt. 
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Figure  3.  -  Comparison  of  the  normalized  brighuiess  of  Mars  (Mariner  9,  Thorpe) 
and  the  value  given  by  equation  (1). 


Figure  4.  -  Comparison  of  the  Minnaert  exponent  for  Mars  (Mariner  9,  Thorpe) 
and  the  value  given  by  equation  (3). 


